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We study the thermodynamics of ultrasmall metallic grains with level spacing <5 compa- 
rable or smaller than the pairing correlation energy, at finite temperatures, T>5. We 
describe a method which allows to find quantum corrections to the effect of classical 
fluctuations. We present results for thermodynamic quantities in ordered grains and for 
the reentrant odd susceptibility in disordered grains. 

1. Outline of the problem 

Superconducting coherence in finite systems has recently received a renewed 
attention since the series of experiments of Ralph Black and Tinkham (RBT).lil 
They succeeded to perform tunneling spectroscopy on individual nanoscale Al grains 
where Ec, the charging energyta, is very large (up to 500 if), enforcing the conser- 
vation of the total number of electrons N. They observed spectroscopic features 
related to superconductivity, which also depended on the parity (N being odd or 
even) of the grains. 

The question of the size limit for a metal particle to have superconducting prop- 
erties was posed since the early days of BCS theory by AndersonE It was argued 
that when the average level spacing 6 ~ l/(J\f(0)V) becomes of the order of the 
BCS gap A, superconductivity should disappear. 

In small metallic grainsEl thermal fluctuations wash out the hallmarks of Cooper 
pair condensation, like the zero resistance and the complete Meissner effect. How- 
ever other features related to superconductivity can be clearly observed. In mi- 
crometres size samples charging effects determine the Coulomb Blockadea and parity 
effectsBQ, but superconducting correlations are still strong (the typical energy scale 
is the BCS A). For ultrasmall grains (S > A, "superconducting dots") at low tem- 
perature, quantum fluctuations due to both large Ec and finite 5 tend to suppress 
superconducting properties. n 

Theoretical analysis by parity dependent mean field theoriesEl found that the 
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breakdown of BCS superconductivity occurs at a value of 8/ A which depends on 
parity. The use of the mean field grand-canonical approach is appropriate when the 
grains are large enough to give small relative fluctuations of the electron number 
(this condition is not strictly met in the experimental!) and it breaks down com- 
pletely in superconducting dots at T = 0. In the canonical ensemble the BCS pair 
amplitude is zero and m order to characterize pairing correlations different quanti- 
ties have to be studiedO, for jitjtance the odd-even staggering between the ground 
state energies at different JVBBEij They studied the BCS Hamiltonian 

n,a 

where {to, n} span a shell of f2 doubly degenerate time reversed (cr = ±) single 
particle energy levels (energy e m and annihilation operator c TOj(T ), and A is the BCS 
coupling constant. The effect of a magnetic field H is accounted for by the Zeeman 
term (jig is the Bohr magneton). 

For equally spaced e m = mS, it was found that even is, ultrasmall grains pairing 
correlations manifest with non-perturbative fluctuationsa and that the low-energy 
physics can be expressed by universal functions of 8/ A (A plays the role of an energy 
scale analogous to the Rondo temperature), allowing a quantitative description of 
the full crossover regime.! ■— ■ 

The model Eq.([l|) was further studied with various methodalil until it was re& 
ognized that it has been exactly solved long ago by Richardson and ShermanJl3 
The ccpiiections of this solution with conformal field theories and off-shell Bcthc 
AnsatzcJ has been recently exploited, and may yield exact expressions for the cor- 
relation functions. 

It is in general very difficult to achieve experimentally unambiguous evidence 
of the presence of strong pairing correlations in ultrasmallr-grains. However the 
spin susceptibility in grains with iV odd Xo{T) was foundo to show a striking 
reentrant behavior (see Fig.l). This qualitative difference with normal metal grains 
is a signature of pairing and persists even in ultrasmall grains. It is due to the 
interplay of pairing which tends to suppress Xo{T) for decreasing T, and parity 
which determines a Curie 1/T contribution which eventually prevails. 

The result was obtained by combinins-a parity-projected path-integral approach 
in the static-path approximation (SPA)Ej with low temperature results obtained 
by massive use of the exact solution. In this paper we describe how to study 
quantum corrections to the SPA with a functional procedure which has been named 
RPA'. The RPA' also allows to study the thermodynamics of weakly disordered 
grains at finite temperature, in an efficient and reliable way. This is an important 
virtue because experiments should be carried out in ensembles of grains, where a 
distribution of shape and volume has to be considered. 

We remark that the Hamiltonian Eq. ([I]) .describes universal properties of small 
grains in the metallic regime (weak disorder if we choose e m as the eigenvalue of 
a random matrb£2] ofj-jthe Gaussian Orthogonal Ensemble (GOE), since we ignore 
here spin-orbit effects.Q A set of mono dispersed grains should be studied along these 
lines. This makes practically impossible the use of the exact solution as in the work 
of Di Lorenzo et al.li-1 Some effects of disorder have been previously studied »sing 
a parity-projected mean field approachta valid for large grains, and at T = OxB 

To clarify the virtues and differences of various approaches we sketch the pro- 
cedure we employ. We start from a path-integral formulation of the problem which 
we consider the grand-canonical ensembles of odd or even grains We express the 
partition function as a path integral over a quantum auxiliary field A(t). The SPA 

consists in retaining only the contribution of all constant A(t) = A", which is 



^ C m, + C m,~ c ri,~ c n,+ ■ (1) 

m,n— 1 
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maximal if = A, the BCS gap. Thus the SPA describes exactly classical fluc- 
tuations. We account for quantum fluctuations by allowing for non constant paths 
A(t). The RPA' consists in retaining Gaussian fluctuations along the static path 

by first fipflgi A j then evaluating Gaussian fluctuations and finally integrating 
over A( s )£3'EiH23 More accurate approximation schemes are also possible. All these 
methods break down at low temperatures, in a more or less spectacular way. In this 
paper we use study the thermodynamics at finite temperature beyond the SPA. We 
present results on the specific heat and on the odd spin susceptibility for equally 
spaced e n . We finally check that the reentrant behavior of \ persists in disordered 
grains. 



2. Functional formulation 

The spin susceptibility and the specific heat of a grain with an even (e) or an odd 
(o) number N of electrons is defined starting from the partition functions Z(T, N), 
which should be evaluated in the canonical ensemble. With the help of a parity 
projection techniqudj the grand partition function of an even or odd system can be 
expressed as 

Z e/o (T, M) = 2 E e " N/T ^ ± ^ N ]Z(T, N) 

N=0 

= \{Z + ±Z_). (2) 

Here, the partition function Z + is the usual grand partition function at temperature 
T and chemical potential /i + = /i. The grand partition function Z_ describes an 
auxiliary ensemble at temperature T and chemical potential /i_ = jj, + itrT] it is 
a formal tool, necessary to include parity effects. In the calculation, the chemical 
potential fj, will be placed between the topmost occupied level and the lowest un- 
occupied level in the even case, while it will be at the singly occupied level in the 
odd case. 

The method we use to perform the calculation is based on the functional formal- 
ism and the Hubbard-Stratonovic (HS) approach, where the two-body interaction 

c m +<y m _c ny -C n ^+ is reduced to a one-body interaction with a time dependent ex- 
ternal field A(r), which has to be averaged over a Gaussian weight. So that the 
problem switches to the treatment of the interaction vertex between the fermionic 
and the HS field. The functional integral over the fermionic degrees of freedom is 
Gaussianr-apd can be easily performed. After this step the grand partition functions 
Z± readsB 

Z ± =M(^Ac,J(dr[j:(TrlnG m ] ± - ^ - " + - 1AJ!]} , (3) 

o 

where, (3 — 1/T and N is the normalization factor of the Gaussian "measure" 
introduced by the HS procedure. The matrix Green function G m> ± is given by 

'-d T - e m + n + n B H A(t) \ 

<W=| U(r-r'). (4) 

A t (t) -d T + e m -fi + fj, B Hl 



It is worth to stress that the {±} ensembles are characterized by different time- 
periodicity of the Green function: in the { — } case the adding term to the chemical 
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potential iirT has been re-absorbed by the gauge transformation on the fermionic 
variables c mi0 -(r) = exp(iirTT)dk a changing, in this way, the periodic boundary 
condition. In particular the fermionic (bosonic) Matsubara periodicity is associated 
to the Z + {ZJ) partition function. 

In order to deal with the non- linear functional integral (^), it is useful to extract 
the static term by the HS field, i.e. A(r) = A(u ) + S n ^o A(w„) exp(— ioj n t) = 
A( s ) + <5A(t). The SPaEI approximation consists in taking into account the only 
static contribution A^, so that the functional integral is therefore reduced to an 
ordinary integral over the class of "straight paths" in the time interval [0,(3]. It 

is worth to notice that A^ s ) is just the average point of the path (the centroid in 
FeynmanJanguage) and for this reason some analogies. with the Feynman variational 
approaches can be found. As in the variational caseo'EEl, the SPA can be improved 
by taking into account small-amplitude (quadratic) fluctuation around the centroid. 
This ifethe so-called RPA', introduced by Wang et a/H3 and refined later by Keiter 
et aLcll in order to study the Anderson model (for an exhaustive review on this 
approximation see also Ref.c3) . In particular it turns out that RPA' corresponds to 
a ring-diagram summation equivalent to the standard random phase approximation 
(RPA), except that only non-zero- frequency transfers at the vertices are retai ned in 
the ring-diagrams, while all zero-frequency transfers are treated exactlyr 1 rr 2 l 
The logarithm term in Eq. (^) can be expanded in the following way 

lnG"! ± = lnG^ 1 + £ i^_( G « ± pf , (5) 

1=1 1 

where 

Performing the trace operation and paying attention to the fact that the zero fre- 
quency are ruled out from the sums, it is easy to find the first order contribution 
vanish. So the static path A( s ) is a local extreme of the integrand. The leading 
correction to the SPA is quadratic in the fluctuation around the centroid. In RPA' 
the partition function reads 

/>oo 

Z±=Af dA^ exp [^± A ( A< s > ) + 8\ ± (A^ )] , (6) 
Jo 

where 



with = ^ + s and £ m = e m — fx, is the standard SPA contributiorJiil't^l. The 
quantum fluctuation correction to the free energy is 

exp [<5V±(A (s) )] = J V 2 [5A]exp{ - 5 2 S±[8A(t)}} , 

where the funcional 8 2 S±[6A] containins products of two Green's functions. The 
path integral can be easily evaluated after the standard sum procedure over (fermionic 



Jcosh *™ + °T \ 

ln2 Unh g " + 2 ^ g J 
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and bosonic) Matsubara frequency has been performed. After some algebra the fol- 
lowing expression is found 



ui n >0 



(8) 



where 

A±(u n )= 2jff m (a; n )taiih 



bm 1 rn 



2T + A' 9mlW " j_ 2^ m (^ + 4^) ' 



(9) 

By means of expressions (^[^) the thermodynamical quantity of the model can 
be numerically evaluated with a CPU time much shorter than what needed for the 
exact solution. This makes possible the analysis of the effect of a statistical distri- 
bution of non interacting levels. In the next section we will show our results for 
even specific heat and for the reentrant susceptibility in the ordered level seating 
case and we compare them with previous results both analytic and exact O This 
comparison, besides, allows to test the reliability of the RPA' in dealing with the 
quantum fluctuations in ultrasmall grains. Finally we will show results for disor- 
dered superconducting dots. 

3. Results 

Results obtained with RPA' are shown in Figs. 1-3. We notice first of all that 
for Xo(T) (Fig.l) the RPA' curves (full lines) move away from the SPA curves (dot- 
dashed) towards the exact result (dashed curve). So the quantum corrections are 
appreciable and give the expected trend. The difference between RPA' and exact 
curves can be mainly attributed to differeiajces between canonical ensemble and 
paruty-projected grand canonical ensemble.Ej The agreement at lower temperature 




Fig. 1. Spin susceptibility for an odd grain vs. temperature in units of 8. Solid line: RPA'; dashed 
line: exact; dot-dashed: SPA 
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Fig. 2. (a) Specific heath for an even grain vs. temperature in units of <5. (b) Rescaled specific 
heat CvS/'yT vs. temperature in units of 8 (see text.). Solid curves are the RPA' results whereas 
dashed curves are the exact results. Thin lines represent the asymptotic behavior at large T in 
the canonical (dashed) and in the grandcanonocal (solid) ensemble. 



is due to the fact that essentially only the classical Curie contribution is important. 
At very low temperature the SPA gives wrong results for the contribution of the 
"condensate" which is by itself small, whereas the RPA' breaks down. For even 
grains Xe(T) vanishes exponentially for T — » 0. 

The specific heat for even grains is shown in Fig. 2a. For large enough coupling an 
anomaly develops. This is not the case for odd grains, as expected, because the odd 
electron reduces the phase space for pairing to be active. At lower temperatures 
the RPA' data seem to reach the data obtained by diagonalization with a good 
quantitative accuracy, at even lower temperatures RPA' breaks down. While is not 
surprising that the Richardson and the parity projected RPA' curves do not match 
each other at high temperature: the former arise from a correct canonical analysis, 
the latter come from the parity projection technique, which relies on the grand- 
canonical treatment. As shown by Denton et alEQ for a normal metal particle, 
the specific heat in the grand-canonical case is lager than in the canonical case, 
asymptotically is ^fce higher, because there are more excitations allowed when the 
electron number is not conserved. In order to point out this behavior, in Fig. 2b 
c v 5/"fT is shown, where 7 = 27r 2 /3 is the Sommcrfeld coefficient of the specific heat 
for an electron gas in the grand-canonical ensemble. 

Experimental investigation on thermodynamic properties of ultrasmall super- 
conducting grains can only be performed in large ensembles rather than in individ- 
ual grains. In the case of the susceptibility it is of course not possible to prepare an 
ensemble of all odd grains. However at temperatures where the Curie contribution 
prevails the susceptibility of even grains is already exponentially small. A similar 
line of reasoning applies for the measurement of the even specific heat anomaly, 
which however seems to be less promising. Two questions can be raised for what 
concerns the physics of ensembles of grains. The first one concerns interaction 
between neighboring grains with possible electron tunneling. However for small 
enough grains large charging energies, as the ones in the RBTlI experiment, prevent 
charge transfer, except for very special situations. Another question is the effect of 
the disorder. To discuss this issue we considered samples of monodispersed grains, 
i.e. grains differ in the shape but not in the volume. This has striking consequences 
in small metallic grains, since shape irregularities over lengths ~ Xp are enough 
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Fig. 3. Effects of a GOE statistical distribution of energy levels on the reentrant behavior of 
Xo(T). Solid lines represent the canonical (A = 0) and RPA' (S/A = 1) results for equally spaced 
noninteracting single-particle levels, while dots are the result for GOE distributed single-particle 
levels 



to affect the spectrum at the Fermi energy. Gorkov and EliashbergEzi incorporated 
the effect of the disordered boundary conditions into the Hamiltonian, where they 
appear as random matrix elements which determine a statistical distributiop-pf the 
energy level. For vanishing magnetic field the proper ensemble is the GOHiZI. The 
thermodynamics of normal disordered grains was studied by Denton et alEB 

We applied the RPA' to this problem. Results are shown in Fig2., where points 
represent ensemble averages over 200 samples of random Hamiltonians belonging 
to the GOE. The spectra were obtained by diagonalizing 500 x 500 matrices with 
random elements and taking the central 100 levels of the resulting semicircular 
distribution, in order to prevent undesired border effects. Results show that Xo(T) 
is clearly reentrant even in disordered samples. 



4. Conclusions 

We have discussed some thermodynamic properties of ultrasmall superconduct- 
ing grains at T > 6 using the RPA' approximation. Comparison with results in the 
SPA, which accounts exactly for classical fluctuations, allows to estimate the effect 
of quantum fluctuations. Comparison with the behavior of the exact solution results 
at T ~ (5 allows to estimate differences between results in the grand canonical and 
in the canonical ensemble. 

The RPA' turns out to be an efficient and reliable method. It can be successfully 
applied to problems for which the exact solution is not useful in practice. An 
example is the study of the effect of a statistical distribution of noninteracting levels 
on the thermodynamics. We considered the odd susceptibility Xo(T) which was 
found to show a characteristic reentrant behavior for grains with a regular spectrum. 
This is a promising quantity for experimental detections of paring correlations in 
ultrasmall grains. We checked that the reentrant behavior is well preserved even in 
samples containing a large number of monodispersed grains, whose spectrum has 
statistical properties described by the GOE. 
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